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Photoassociation spectroscopy experiments on ''"Ca atoms close to the dissociation limit 
4s4s^So— 4s4p^Pi are presented. The vibronic spectrum was measured for detunings of the photoas- 
sociation laser ranging from 0.6 GHz to 68 GHz with respect to the atomic resonance. In contrast 
to previous measurements the rotational splitting of the vibrational lines was fully resolved. Full 
quantum mechanical numerical simulations of the photoassociation spectrum were performed which 
allowed us to put constraints on the possible range of the calcium scattering length to between 50 ao 
and 300 ao. 

PACS numbers: 34.50.Rk, 32.70.Cs, 32.80.Pj, 34.10.+X 

I. INTRODUCTION 

In photoassociation spectroscopy experiments within cold ensembles of atoms a bound molecule is formed under 
the influence of light. Such measurements can yield information about many atomic properties such as excited states 
lifetimes, ground state scattering lengths and long range potentials. Most of the experiments are performed with 
alkaline elements. However, due to their atomic fine and hyperfine structure theoretical calculations are in general 
very complicated to perform for these elements. In contrast, the most abundant isotopes of the alkaline earth metals 
show a non-degenerate ground state with no hyperfine splitting. This makes the comparison between experiment 
and theory easier than in the case of alkali metals. In Fig. |^a a simplified molecular level scheme is shown for 
the case of Ca2. The photoassociation connects the ground state asymptote to the molecular state ^E^. A/ is 
the frequency difference between the photoassociation laser and the asymptote ^Sq + 4p^Pi. Our first observations 
on photoassociation of Ca ^ were performed with cold ensembles around 3 mK and with short trap lifetimes of 
about 50 ms only, which in total leads to a poor signal-to-noise ratio for the trap loss by photoassociation. Only for 
the deepest recorded vibrational level at the excited asymptote we were able to observe rotational structure on the 
vibrational transition. The line profile and line position were interpreted by the help of assumptions, like an averaged 
fixed frequency shift of all transitions by 50 MHz due to 3 mK samples. 

Here we will report of a more complete series of measurements under improved conditions, doubling the spectral in- 
terval of measured photoassociation lines, giving clear rotational structure and performing a full line profile simulation 
to get the proper shift of the line center due to thermal averaging and the intensity of the association laser. 

II. EXPERIMENT 

The calcium atoms are cooled and trapped in a magnetooptical trap (MOT) using the strongly allowed transition 
at 422.79 nm (F = 2.16 • 10^ s^^, see the level scheme in Fig. 0b). An overview of the experimental setup is shown 
in Fig. H The 423 nm radiation is generated by frequency doubling about 160 mW of a master-slave diode laser 
combination with a KNbOs crystal inside an enhancement cavity. The output power at 423 nm is 36 mW with a 
linewidth of 300 kHz. A separate thermal Ca beam is used to keep the absolute frequency in resonance with the 
atomic transition. The MOT consists of three mutually orthogonal retrorefiected beams with diameter of 10 mm and 
a power of 5 mW in each beam. The frequency of the trapping laser is red detuned by 29 MHz with respect to the 
resonance. The atoms are directly captured from the low velocity tail of the Boltzmann distribution of a thermal 
atomic beam. The magnetic field gradient generated by a pair of anti Helmholtz coils is 9 • 10~^ T/cm. A repump laser 
at 672 nm is used to prevent trap losses via the ^Di state, thereby increasing the density of the atomic ensemble by a 
factor of four. Absorption imaging of the cloud on the cooling transition was used to determine the radius r and the 
density no of the trapped ensemble leading to typical values of r = 0.3 mm and no = 2 • 10^" cm^^. The temperature 
of the ensemble was deduced independently from the ballistic expansion of the cloud and from the Doppler width of 
the excitation spectrum on the narrow ^So — '^Pi intercombination line (F = 2000 s~^), yielding 2 mK < T < 3 mK. 
With the repump laser turned on the lifetime of the trap is around 300 ms. 

The light of the photoassociation laser is generated by a Stilben 3 dye laser pumped by a 3 W Ar+ UV laser. 
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The output power is 100 mW at 423 nm with a linewidth of 1 MHz. The laser is frequency stabilized to a tunable 
high finesse cavity in vacuum. A tunable hehum neon laser at 633 nm is stabihzed to the same cavity and its 
beat frequency with an iodine stabilized hehum neon laser is measured. By keeping the beat frequency constant 
the frequency of the photoassociation laser can be held fixed with an accuracy of ±2 MHz. Furthermore, with this 
scheme it is possible to scan the frequency of the photoassociation laser by approx. 1.5 GHz with the same accuracy. 
The detuning A/ — fca. — /Laser of the photoassociation laser with respect to the atomic resonance ^So — ^Pi was 
determined by measuring its beat frequency with the frequency stabilized trapping laser for detunings up to 20 GHz. 
For larger detunings the frequency of the photoassociation laser is determined by means of a wavemeter which allows 
absolute frequency measurements with an accuracy better than 0.2 MHz. With the transition frequency of the atomic 
resonance determined the same way the detuning can be calculated. 

The rate of photoassociation was measured by recording the increase of trap losses due to photoassociation which 
manifest themselves in a reduced trap fluorescence. The state of the photoassociated Ca dimer may undergo a state 
change to a molecular state dissociating into atoms which are no longer trapped (state changing collision), or the 
dissociative decay of the dimer may lead to a gain in kinetic energy high enough to allow the atoms of the pair to escape 
from the trap (radiative escape). The photoassociation laser was switched on and off by means of an acousto optical 
modulator (AOM) for time intervals of 500 ms which must be considerably larger than the trap lifetime of 300 ms 
to allow the trap to reach a new steady state. The intensity of the photoassociation laser was around 400 mW/cm^. 
Before and after each interaction period with the photoassociation laser the 423 nm trap fluorescence was detected 
by an avalanche photodiode. The coUisional loss fraction is the ratio of the reduced signal due to photoassociation 
and the signal obtained without photoassociation. The trapping laser itself stayed on during the whole time. Due 
to intensity fluctuations of the trap fluorescence of 1 % to 2 % averaging is done over at least 30 s to increase the 
signal-to-noise ratio. The coUisional loss fraction is directly proportional to the atomic density of the ensemble and, 
as long as saturation effects can be neglected, proportional to the intensity of the photoassociation laser. 

We measured the coUisional trap loss for detunings of the photoassociation laser ranging from 0.6 GHz to 68 GHz 
below the atomic ^So — ^Pi transition. If counting the vibrational levels v' from the asymptote starting at w' = 1 
this frequency interval corresponds to 34 < w' < 72. To smaller v' we were limited by the resonant interaction of 
the photoassociation beam with single atoms, pushing them out of the trap by radiation pressure. For higher v' the 
measurements were limited by decreasing signal strength. 



III. THEORETICAL DESCRIPTION 



The theoretical description of photoassociation spectra consists of two parts: the precise determination of the 
bound states close to the dissociation limit and the modelling of bound-free transition together with its detection by 
trap loss. To describe the asymptotic level structure for a diatomic molecule like Ca2 we use the accumulated phase 
approximation which was described by VerHaar et al ||^ . The potential is given by the long range branch according to 
the atomic states at the asymptote. For Ca2 at the excited asymptote ^So-f ^Pi the long range behaviour is dominated 
by the resonant dipole-dipole interaction and the centrifugal potential term for the rotational state J which should 
have the correct behaviour for dissociating into an atom pair with non-vanishing atomic angular momentum; for the 
molecular state it becomes: 

U(R)-n Ca ft^[J(J+l)+2] 

UiR)~D-- + (1) 

Here D is the asymptotic atom pair energy and C3 is the dipole-dipole interaction constant which can be calculated 
from the transition moment of the "'^Sq — > "'^Pi transition with frequency uja or from the natural linewidth 7^ (all 
linewidths in this paper are in energy units) of this transition, giving C3 — 3c^7a/2aja^ for the molecular state -^T,^, 
and fi is the reduced mass of the molecule. Because of applying only the potential branch at large internuclear 
separations we replace the left boundary condition (small R) for a bound state by a fixed phase which accumulates 
from small internuclear separation up to a large R, say Rq, which is still small enough compared to the classical outer 
turning point of the desired bound level. If the binding energy of asymptotic states is small compared to the kinetic 
energy for the vibrational motion around the minimum of the potential the accumulated phase up to Rq is a very 
weakly varying function for all such asymptotic levels. Thus a fit of observed levels measured with respect to D can 
be performed with the above potential branch and a simple power expansion of the phase with the binding energy as 
expansion coefficient. It is numerically easier to choose a weak variation of Rq for a constant phase, where in our case 
we programmed a maximum of the wave function at the slowly varying Rq. 

For modelling the photoassociation profile we apply the theory developed by Bohn and Julienne |^ . The coupling 
of the continuum at the ground state asymptote ^Sq-I-^So by the laser field to the excited asymptote ^Sg-l-^Pi leads 
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to an excitation rate F and an energy shift Ei of the transition. The above paper shows that the trap loss probabihty 
by photoassociation for a coUision pair with a kinetic energy e can be written in the following form: 

|o ils + 7r)r 

I o'l [e-(A-£;i)]^ + ( ^°+^-+^"+^ )^ 

Sol is the scattering matrix element of the process of two colliding atoms on the ground state via the bound intermediate 
state by the light field to the loss state /. 7s and •jr are the rates for state changing collisions and the radiative escape, 
respectively. 7ts is the molecular spontaneous decay back to trapped states. A represents the detuning between the 
excited molecular bound state with energy Eb{v, J) and the laser frequency with respect to the ground state asymptote: 
A ~ Ef, — hujL (see Fig. ||). A must not be confused with A/ which gives the detuning of the photoassociation laser 
(see Fig. |l|a). In the theoretical description we will count the vibrational levels v from the bottom of the potential 
well. With vd &s the non- integer vibrational 'quantum number' derived from the vibrational phase of the hypothetical 
state exactly at the asymptote one gets the relation v = vd — {vo modulo 1) — u' + 1. F and Ei can be evaluated 
in the reflection approximation, where the scattering wave function of the ground state potential for the scattering 
energy e is only needed at the Condon point of the transition, the above mentioned paper gives the results: 

r-MK,.)^— ^ (3) 

p ,^ M k{Rc)9{Rc) 
E,=.iVob)^ • 

Here the derivative of Eh with respect to v gives the vibrational spacing around the considered bound state, and fo 
and go are the regular and irregular solution of the scattering function for kinetic energy e. Rc is the Condon point 
of the potentials interacting by the laser field at the detuning A. Dc is the potential gradient of the excited state 
at the Condon point, the contribution from the ground state potential plays no rule at large internuclear separation 
considered in the present case, because this behaves like 1/R^ compared to 1/R^ for the excited state. Vb& is the laser 
interaction which is calculated with the help of the laser intensity /, the atomic natural width 7^ and the molecular 
enhancement factor /moi which is determined by the molecular symmetry and of the transition and is to 
good approximation here 2 (see for example the dipole formula in [Q): 

(Vob) = 5 ^ 

LOa is the atomic transition frequency, here of the ^So to ^Pi transition. The photoassociation rate K{A.,T) at 
the detuning A from a thermal ensemble at temperature T is then obtained by thermal averaging of the trap loss 
probability and summing over all contributing partial waves € to the bound state J: 



and 



1 r°° 



Sol is shown with all dependencies given in the equations above and Qt is the partition function for T for the atom 
pairs with reduced mass fj,: 



Qt = 



/27r/ifcr\ 2 



For the simulations of photoassociation profiles, especially for the comparison of relative intensities between different 
vibrational states, one needs a reliable estimate of 7^ and 7^ at least the relative function with vibrational states. 
The losses can be either state changing collisions or radiative escape processes, for both we need an inward motion 
of the photoassociated state. Thus an estimate of the relative probability of the loss will be the vibrational period of 
the state if one assumes that each crossing of the inner region of the potential leads to an equal amount of loss. The 
vibrational period can be estimated from the LeRoy-Bernstein formula for bound states close to the asymptote: 
Eh = D — X(vd — vY where X is a constant depending on C„ of the long range behaviour and the reduced mass of 
the molecule. In the case of the dipole-dipole asymptote C3 is involved and i becomes 6. With this approximation the 
loss rate 7s + 7^ becomes proportional to 6X{vd — v)^; this value was introduced into the simulations of the profiles. 
For calculating the eigenenergies with the accumulated phase approximation and the scattering wave function the 
conventional Numerov method was applied with sufficiently small integration steps, which was always checked by 
numerical convergence. 
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IV. EXPERIMENTAL RESULTS 



Figure ^ shows the trap loss spectrum for the vibrational levels v' — 69 and 72 of the state at a detuning 
A/ of the photoassociation laser around 52 GHz and 68 GHz from the atomic resonance. Three rotational lines 
corresponding to J = 1, 3 and 5 are observed. The sample temperature is about 2 mK, at which the collisional partial 
waves with € > 6 have too low probability for leading to an observable signal of level J = 7 of the excited state. 
The solid curve shows a simulated spectrum. It is obtained by starting with a potential that is known, e.g., from our 
earlier spectroscopic measurements 1 1[ , from which the eigenenergies of the rovibrational levels are calculated with the 
accumulated phase method. These do not necessarily coincide with the positions of the vibrational levels as measured 
in the experiments due to the finite temperature, AC Stark shifts and the potential itself. The position of rotational 
levels with J values of 3 and 5 are shifted up to 120 MHz by these effects. The shift between the observed spectrum and 
the simulated one is used in an iteration procedure to modify the potential described by C3 and the power expansion 
of Rq. After two to three iterations convergence is achieved. The procedure gives a C3 = 0.5217(45) ■ 10"^ nm^ ■ cm~^ 
resulting in an atomic decay rate of 7a/^ = 2.150(19) • 10^ s^^. This value is smaller than that reported in |jl| by three 
times the error limit. This can be explained by the fact that in ^ shifts of the line positions due to temperature have 
only been accounted for by a global shift of all frequencies of 50 MHz. From our new calculations presented here this 
is not justified and thus leads to a different value. It was also checked that the retardation effect of the dipole-dipole 
interaction must be included in the fit. Especially, the levels close to the asymptote ^So+^Pi for v' < 40 otherwise 
show non-negligible deviations. 

It turns out that the relative height of the components of the rotational substructure depends significantly on the 
value of the scattering length for the ground state potential. For varying the scattering length in test calculations, 
one may change the potential slightly at the repulsive branch at small R or around the minimum of the potential. 
The result are slight shifts of the last bound levels not only for £ — but also ioi £ ^ 0. Here for £ = 2 the position 
of the corresponding shape resonance moves to the asymptote for decreasing scattering length, which results in the 
variation of the intensity ratio of the peaks for J — 1 and 3. Fig. ^ shows the numerical simulation of the vibrational 
level v' = 69 for three different scattering lengths of 52.4 cq, 85.2 and 324 gq, respectively, with the Bohr radius 
oo = 0.0529 . . . nm. For clarity the curves are shifted in the vertical direction. The ratio of the J = 1 and J = 3 
components and the deepness of the dip between the two components of all measured vibrational lines with resolved 
rotational substructure {v' > 62) are best reproduced by the potential with a scattering length of 85.2 oq. The quite 
different form of the two other curves (52.4 ao and 324 ao) suggests that the scattering length lies somewhere between 
50 Oq and 300 oq. This is the same order of magnitude of the scattering length as determined by other spectroscopic 
methods 

In our previous photoassociation experiments on '*'^Ca it was not possible to determine the relative intensities of 
the vibrational lines. This was mainly due to the short trap lifetime of 50 ms which leads to maximum collisional trap 
loss fractions of 1 %. In the experiment reported here, the collisional loss fraction at the center of each vibrational 
line was measured (Fig. For vibrational lines where the rotational structure is not resolved the collisional loss 
fraction at the expected maximum of the resonance line was determined (full circles). For large v' the rotational 
structure is sufficiently resolved and the collisional loss fraction was measured at the position of the J = 1 (full 
triangles) and the J = 3 rotational lines (open triangles). The collisional loss fraction with the photoassociation laser 
tuned to the gap between two vibrational lines is shown by open circles. For v' < 47 the vibrational lines begin to 
overlap and the collisional loss fraction does not drop to zero between two lines. Before and after each measurement 
session of approximately ten hours the atomic density was determined by taking absorption images on the cooling 
transition. Due to mechanical drifts in the setup the density decreased over this period of time by up to 20 %. A 
linear variation of the atomic density with time was taken into account to normalize the collisional loss fraction to a 
density of 2 ■ 10^° cm^^. In addition the change in the intensity of the photoassociation laser that showed fluctuations 
of approx. 10 % during a scan was corrected. From the measured variation of collisional trap losses with the intensity 
of the photoassociation laser (Fig. ^ the fluctuations were corrected to an intensity of 280 mW/cm^. 

The curves in Fig. ^ depict the intensity variation of the trap loss as a function of the vibrational quantum number 
v' as expected from theory where the functional form of the trap loss in its vibrational dependence is simplified being 
proportional to {vd — v)^ from the local vibrational spacing. This approach does not distinguish the contribution by 
state changing collisions from those by radiative escape. Thus at least a two-parameter ansatz would be necessary 
for which a functional model does not exist in the literature for describing the vibrational dependence. We were 
experimenting with different exponents for (vd ^ instead of 5 and with adding a constant to the loss rate as a 
second parameter which could model the different weights of the two mechanisms. Such an approach gives a nicer 
looking graph for the comparison between experiment and theory but not a convincing deeper physical insight. Thus 
we suppressed it, but we would prefer to indicate clearly with Fig. ^ the status of our present understanding. 
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a) b) 




FIG. 1: Excerpt from the Ca molecular (a) and atomic (b) level scheme. The transition around A = 423 nm is used for cooling 
and photoassociation. A repump laser (A = 672 nm) prevents decay into metastable triplet levels via the ^ D2 state. The long 
range part of the ^T,^ molecular potential is probed in the experiment. The potential energies as a function of the internuclear 
distance 7? are indicated schematically. A/ is the frequency difference between the photoassociation laser and the atomic 
resonance. 



V. CONCLUSION 



In conclusion, photoassociation measurements on cold calcium atoms for an extended range of detunings of the 
photoassociation laser ranging from 0.6 GHz to 68 GHz below the asymptote were performed and compared with 
full line profile simulations. We find good agreement between the observed and calculated spectra. The observed 
intensity ratio of J = 1 and J = 3 resonances determines the range of the ground state scattering length between 
50 ttQ and 300 ag taking the long range behaviour as determined by spectroscopic experiments or by theoretical 
predictions . The measured relative intensities of the vibrational transitions are not well explained. New studies of 
loss mechanisms are needed or better electronic structure calculations on Ga2 would be desirable to improve on the 
model calculations which were published by Machholm et al. Q. 

Up to now we only observed the creation of calcium dimers in the ^Y.^ state. Besides this, the formation of atom 
pairs in the ^Hg state is allowed due to retardation. Those vibrational levels should show up as very sharp features 
(a few MHz width) in the spectra It may be possible to get signatures from these states in the experiment when 
switching to ultracold calcium atoms with a temperature of a few microkelvin. 
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This work was supported by the Deutsche Forschungsgemeinschaft. 
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FIG. 2: The setup used in the experiments. Calcium atoms are captured from a thermal beam and cooled in a magnetooptical 
trap with three mutually orthogonal, retroreflected beams in cr+ — a- configuration (vertical beams not shown). The 423 nm 
trap light is generated by a frequency doubled master-slave diode laser combination. A dye laser produced the light used for 
photoassociation. For details see text. 
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FIG. 3: Schematic potential diagram illustrating the notation used in the text for describing cold collision processes. 
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FIG. 4: Vibrational lines v' = 69 and 72 with three resolved rotational lines corresponding to J = 1, 3 and 5. The solid curve 
represents the prediction of a full quantum mechanical calculation assuming a scattering length of 85.2 oo. 
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FIG. 5: Simulated spectra for the vibrational level v' = 69 for 2 mK and 400 mW/cm^ photoassociation laser intensity for 
three different values of the scattering length. The curves are shifted in the vertical direction for clarity. 
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FIG. 6: Measured coUisional trap loss as a function of the vibrational quantum number v' . For vibrational lines where the 
rotational substructure overlaps the loss was determined at the expected maximum of the resonance line (full circles). For lines 
with sufficiently resolved rotational structure the loss was measured for J = 1 (full triangles) and J = 3 (open triangles). The 
curves indicate theoretical predictions with simplified assumptions about decay rates (see text). 
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FIG. 7: Dependence of the coUisional loss fraction at the peak position of line v' = 39 (open circles) and line v' = 49 (full 
circles) on the photoassociation laser intensity. At the highest intensities used in the experiments saturation effects begin to 
show up. 



